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Abstract 

The leading, planar diagrams of the 1/N C expansion and the usual string description sug- 
gest that quarks propagate on the boundary of a two-dimensional world surface. We restrict 
the quarks to the boundary of the world surface by giving them infinitely large mass on 
the interior of the surface and zero mass on its boundary and show that in this picture the 
QCD #-vacua can be represented by the self-intersection number (or equivalently by the 
first Chern number of the normal bundle) of the surface. 
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1. Introduction 



In recent years, an effort has been made to establish the relationship between QCD 
and strings in order to better understand QCD. Indications that strings might provide an 
effective description of QCD are seen in the lattice formulation and the 1/N C expansion 
[1,2]. In the limit, N c — > oo, only planar diagrams, which have exactly the topology of 
the quantized string with quarks at its ends, contribute. The two-dimensional structure 
obtained by attaching surfaces to each loop in the planar diagram resembles the world sheet 
E of an open string. If QCD behaves like a string theory, we should find string theoretic 
representations of QCD characteristics. One of the qualitatively important characteristics 
of QCD is the existence of 6>-vacua. The string-candidate associated with the 6>-term of 
QCD, j^trFF, was suggested in [3,4] to be the self-intersection number of the surface, 
/(E) [5]. In [4] Mazur and Nair introduced fermions in the theory and projected the 
four-dimensional action onto E. This projection determines the coupling of the fermions 
to the extrinsic and intrinsic parameters of E. They found that the effect of 6>-vacua, 
represented in the QCD functional integral by exp(i9Q) where Q is the instanton number, 
can be represented in the string description by exp(«6>Ci(E)) where Ci(E) is the first 
Chern number of the normal bundle to E. However in the large iV c -limit description 
or the usual string description the quarks are constrained to the boundary of the two- 
dimensional surface, and in [4] the quarks live on the entire surface. Therefore a more 
accurate description of QCD would require a modified model. In this work we restrict the 
quarks to the boundary by giving them infinitely large mass on the interior of the surface 
and zero mass on the boundary. We calculate the axial anomaly for these quarks with a 
coordinate-dependent mass and obtain a similar relation between the first Chern number 
of the surface and the total change of the axial charge, AQ 5 . 

2. Short review 

Following Mazur and Nair [4], we begin with the four-dimensional free fermion La- 
grangian 
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C = ^ 1 a {d a --wf[ la ^p])^ (1) 

where we have pulled back quantities from E to M and is the spin connection. M 
in general is a Riemann surface with real coordinates £ a (a = 1,2). The image of M in 
1Z 4 is the surface E. More than one point in M can be mapped to the same point in 1Z A 
corresponding to self-intersections of E. Thus the mapping : M — > 1Z 4 is, in general, 

only an immersion, can be expressed in terms of the tangent and normal vectors of 

the immersion X^(Z) to give the Lagrangian for quarks on M. The tangent vectors are 
t% = dX^/d^ a and the normal vectors are orthogonal to the tangents at the image points 
on E. The immersion of M in 1Z 4 defines the tangent and normal bundles on M. A a AB is 
the SO (2) connection for the normal bundle. The extrinsic curvature K^ b is given by the 
equation 

d a d b x» = r c ab d c x» + K*n A 

where T^ b is the usual Christoffel symbol for the induced metric, g ab = dX^/d^ dX^/d^ b . 
The frame fields on M, e", are given by e a a eab = 9ab- The target space is 1Z 4 7 thus the zero- 
curvature spin connection is 

where E a ^ decomposes as follows 

E afl = (n A ^,t^e a& ). 
Hence the Lagrangian for fermions on M may be written as 

C = *l a (d a - ^ 6 [/, 7c ] + ^ aAB [ 7 ^, 7 B ] - \K*[ lb , lA })*. 

Note that this Lagrangian only contains geometric information. The strings description 
of QCD should incorporate confinement, thus to obtain a purely geometrical coupling for 
color singlet combinations there should be an extra factor of 1/N C in the couplings of K^ h 
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and A a AB to quarks. This factor arises naturally in the large iV c -limit of QCD since the 
quark contribution is smaller by a factor of 1/N C relative to the gluon contribution. Thus 
the Lagrangian for quarks on M can be taken to be 



C = q^ a D a q t 

where D a is given by 

Da = d a - ^ b [ 7 6 , 7c ] + ^- A^fy 4 , 7 s ] - ^-^[76, 1a\- (2) 

The standard heat kernel method can be used to obtain the axial anomaly associated with 
the covariant derivative D a . The functional trace of 75, the quantity equal to the instanton 
number in four-dimensional QCD, was found to be [4] 

^75 = tJ-^i(S) (3) 
where Ci(E), the first Chern number, is related to the self intersection number /(£) by 

id(E) = 1(E) = --L J d 2 ^g%V a t^V b t^, 
and has the simple form 

Ci(E) = i- f trF. (4) 

F is the field strength tensor for the SO(2) normal bundle connection and t^ v = ^=t^. 

A complete analogy with QCD is lacking here because in this picture the quarks live on 
the two-dimensional surface, rather then propagating only on the boundary as is required 
by the leading order diagrams of the 1/N C expansion and by the usual string description. 
To provide a better analogy we instead choose a regularization scheme which allows us to 
compute the anomaly in two dimensions, while restricting the quarks to the boundary of 
the surface. 



3.Regularization scheme 



We give the quarks zero mass on the boundary and infinitely large mass inside the 
surface. We shall Calculate the propagator for these coordinate-dependent mass quarks 
and prove that they are indeed restricted to the boundary. For simplicity we take space 
to be flat with the spatial direction a varying between — / and /. m(a), the quark mass, is 
chosen to have the simple form: 

{0, — I < a < —a region I; 
mo, —a < a < a region II ; (5) 
0, a < a < I region III . 

We are interested in the limit mo — > oo. We shall assume that the gauge field of interest, 

viz. {F a b)AB, the field-strength for A a AB, is slowly varying on the scale of (/ — a). 

4 .Calculating the anomaly for a coordinate- dependent mass 

In order to derive an equation similar to equation (3) we calculate the quantity AQ 5 
given by 

AQ 5 = J d^d 2 i (6) 
where J 5m is the neutral gauge invariant axial current, J 5m = vE^^ 5 !]/, and 

=A + 2im(x)y(x)<y 5 -$(x) (7) 

where A is the anomalous term. There are different ways to calculate the axial anomaly. 
The perturbative method is best suited to our choice of mass. The two-dimensional axial 
anomaly stems from the two-point function 

Tl = (T(JM4(x))). 

For a particle with constant mass, m, requiring J^(x) to be gauge invariant results in an 
anomalous term, p^- 

5 



P2 = —iTF^ix). (8) 

Since we want to calculate q v T^ v in a spatially confined space we choose a cutoff regular- 
ization. We are also interested in the limit m — > oo, thus we choose the cutoff A and the 
fixed mass mo such that 

A < m . (9) 

We next compute the propagator for m = m(a). In view of equation (9) we search for the 
solution to the Dirac equation: 

■ d -m(a))^(a,t) =0 

with \E\ < A < mo which satisfy the boundary conditions required for a self-adjoint 
hamiltonian, viz: 

* I *(-07°7 1 *i(-0 = 

(10) 

*ih(07V*iii(0 = o 

and the matching conditions at a = ±a, where is the solution in region i. Separation 
of the time-variable requires a two-phase ansatz: 

■qj = e - lEt ¥ + \a) + e lEt ¥-\a). 
The static functions \l/ ± (a) satisfy 

±E^ ± \(i) = (ap + /3m)^ (±) (a). 

Our representation for the Dirac matrices is a = a 2 , (3 = a 1 . Implementing the matching 
conditions and requiring a normalizable solution when mo — > oo, while fixing the wave 
function in region III (or I) gives a solution of exponentially decaying wave functions in 
regions I (or III) and II, but an oscillating wave function in region III (or I) . The zero flux 
condition at x = I gives quantization of energy 
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2k n (l - a) +cp = 27m {n E Z : \k n \ < A} (11) 
where k 2 = E 2 and sin (p 



m 

For large mo, sin </? gives 



2mo 

where 



fcn = 7i r— -iV<n<iV (12) 

(I - a) + 



\N\ < -(— + (/ -a)). 
1 1 - tt 2m V " 



Fixing \1>i instead of yields 



2/c n (/-a) + (^ = 27r(n + l/2) {n E Z : \k n \ < A}. (13) 

The other boundary condition is automatically fulfilled in the large mo limit. The propa- 
gator is constructed from \E , ^ ± ) (x) 
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S(x,x') =-i6(t'-t) e- tE " (t '-^¥+\x')¥+\x) 

+ie(t-t') e iEn ^- t ^ir\x , )^ir\x). 



n=-N 

For k n satisfying (12) the propagator in regions I and II decays exponentially. Note that 
in region III the propagator describes a Dirac particle in an infinite square potential well 
of width (I - a ) + 5^. 

The propagator of a massless fermion in an infinite square potential well of width 
(/ — a), denoted by S°(x',x), is given by equation (14) with 

^(+)/ \ _ 1 ( sinfc n (x' - I) cosk n (x - I) smk n {x' - I) smk n {x - I) 

n n l-a\ cosk n (x' - l)cosk n (x - I) |p- cosk n (x' - I) smk n (x - I) 

(15) 

where 
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k n = -^- N°<-(l-a) -N°<n<N°. (16) 
I — a TV 

^n(x')^n(x) is obtained from (15) by replacing E n with —E n . S°(x',x) does not have 
the usual x' — x dependence. This is not surprising since in confining the quarks we have 
lost translational invariance. It is straightforward to check that S°(x',x) is indeed the 
propagator. The propagator corresponding to a well of width (/ — a) + differs from 
S°(x',x) by terms of order — which appear in the normalization of the wave function, 



1^41 2 _ n \ 

11 4(/-a) V (l-a)2m ; 

and in the argument of the trigonometric functions 



n 2(1 -a) { (l-a)2m ) ' 
(Though the interval in equation (12) appears different than that in equation (16), we are 
free to choose N = N .) When computing T^ u in region III, no additional factors of m,Q 
appear. This insures that the mo — > oo limit gives the same expression for T^ u as in the 
massless case. Thus we find for region III: 



A = p 2 . (17) 

A similar argument for k n = ^["1^ gives the usual two-dimensional anomaly in 
region I. In region II, using the cutoff regularization, the propagator between a and a' 
decays exponentially. Therefore all diagrams in that region vanish in the large mo limit. 
Also all matrix elements of J 5 = mo 1 ^ 5 ^ vanish in this regularization scheme since the 
factor of mo cannot compensate for the exponential decay. This gives in region II, 

= 0, (18) 

which is expected. A simple way to see this is to recall that for a massive fermion in two 
dimensions we have 
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<9 M J 5M (:r) = p 2 + 2zm^(x)7 5 ^(x). 

When m — > oo the matrix elements of \j>7 5 \I/ behave like ^, resulting in a non- vanishing 
remainder for 2im#7 5 \E'. In fact this remainder is exactly — p 2 - Thus we get for a heavy 
fermion 

lim d.J^ix) = 

which is the same result as (18). 

(To understand this better, recall [6] that another way to assure gauge invariance is by 
introducing a Pauli-Villars regulator field ip with mass M and a regulated axial current. 
Then d^J 5 ^(x) = 2im^(x)^ 5 ^(x) - 2iM^(x)^ 5 i)(x). When M -> oo the regulator field 
contribution gives the axial anomaly.) 
To summarize 

{P2, -I <cr < -a; 
0, -a < a < a; 
P2, a < a <l. 

5. Conclusion 

The total change of the axial charge, AQ 5 , constructed from the axial current is 

AQ 5 = fdt [ dap 2 (20) 

iV C J Jl + III 

where the same argument that led to the covariant derivative of equation (3) results in a 
coefficient of 1/N C . It is this quantity that we want to relate to the first Chern number or 
to the self-intersection number of the surface E. (In [4] a relation to the functional trace 
of 7 5 emerges naturally from the functional integral formalism.) In equation (4) Ci(E) is 
expressed as an integral of the field strength tensor, F^ab- Since the self-intersection num- 
ber is a topologically invariant quantity of the surface E, given a connection configuration 
A a AB, we can smoothly deform it to A' aAB , such that 
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(19) 




We can choose this deformation to be the one for which the field strength tensor F ab AB 
vanishes in region II, i.e. where m(a) 7^ 0. This insures that the normal bundle is trivial 
in region II hence only regions I and III contribute to the self-intersection number. Thus 
we can write 



The integrand is exactly p2- In both equations (20) and (21), the area of integration, which 
is the 'dynamic' area, is limited to a strip along the boundary where the quarks live. This 
allows us to relate the QCD topological charge to the first Chern number in the string 
description 



in the same manner as in ref. [4] , but for a more accurate picture in which the quarks are 
restricted to the boundary. 

Note that if we deform A a AB such that a self-intersection occurs in region II, the gauge 
fields are no longer varying slowly on the scale of (I — a). This must be taken into account 
when deriving the anomaly (e.g. requiring a more careful treatment of the boundaries 
between the regions). Presumably, a non- vanishing divergence of the axial current in 
region II will result and equation (22) will still be valid. 
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